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Abstract 



The Hamilton-Jacobi method is generalized in classical and relativistic 
mechanics. The implications in quantum mechanics are discussed in the 
case of Klein-Gordon equation. We find that the wave functions of Klein- 
Gordon theory can be soundly interpreted as describing the motion of an 
ensemble of particles that move under the action of the electromagnetic 
field alone, without quantum potentials, hidden uninterpreted variables, 
or zero point fields, 
pacs 03.20.+i, 03.20.+p, 03.65.-w 

1 Introduction 

In previous papers (pQ, [2]) we have considered the possibility of reinterpreting 
Klein-Gordon and Dirac's equations, as describing ensembles of particles moving 
under the action of the electromagnetic field, where the number of particles is 
not locally conserved. In order to do that for the Klein-Gordon equation, on 
the basis of some general electrodynamic considerations, we proved that it was 
sound to assume that the field of kinetic four-momentum 



where S is the phase of the wave function; is the electrodynamic four- 
potential; and $ is a function of space-time coordinates, whose nature we left 




is given by 



(1) 



completely undetermined. To deal with Dirac's equation we used a similar 
approach. 

In this paper we determine the nature of the function $. Actually, we prove 
that the principle of relativity requires that we consider a general four-vector 
ujfj,, instead of the gradient of a potential. For this we start with an analysis of 
Helmholtz's theorem on the representation of vectors fields and its extensions 
to n dimensional spaces, in particular to space time|3]. In the first section we 
derive Helmholtz theorem: 

/ = V0 + VxA, (2) 

following a variational approach. This permits us to prove that if the potentials 
are chosen in a manner that (V</> — /) • n = at the boundary of the region 
under consideration, then V</> is the best approximation to /, in the sense of 
the quadratic norm — which we think is an original result. 

Further we apply those conclusions to perform an analysis of the physical 
meaning of Hamilton- Jacobi theory. Taking it out of the mathematical realm 
of canonical transformations, allows us to find its necessary generalizations in 
classical and relativistic mechanics, laying thus a foundation for a complete 
explanation of the field ui^. In the second section we work out a complete 
example of this extension of Hamilton- Jacobi theory in a very simple classical 
situation, just to show how it works, as well as to prove that there are very simple 
and common dynamical scenarios that are being disregarded by the present 
theory and, therefore, that our extension is not only a theoretical speculation, 
but a necessary step to extend the scope and strength of a fundamental part of 
classical mechanics. 

We consider a flux that corresponds to an infinite ensemble of particles with 
the initial velocities: 

v(f, 0) = uj x r, (3) 

where a; is a constant vector. We prove that this flux is invertible, in the sense 
that we can obtain the initial position of the particle that passes through the 
point x at time t. However this flux cannot be derived from a potential in the 
form 

v= , 4 

m 

because 

V x v ^ . 

The generalization of Hamilton- Jacobi theory is even more necessary in the 
case of special relativity because — as we will prove — any field of four-velocities 
that describes the motion of an ensemble of particles has a vorticity, at least 
it describes a flux of free particles. Furthermore: all we have is to consider a 
relativistic ensemble of free particles, where the initial momentum is defined as 

p(f, 0) = lu x r, (5) 
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for a particular observer, to have a system whose evolution cannot be described 
by the Hamilton- Jacobi theory as it is presented in pp. 24-29]. In some way 
we follow Bohr when we state that if there is a reformulation of Hamilton-Jacobi 
theory that widens its range of applicability it should be incorporated into the 
body of mechanical knowledge, no matter if it is or not compatible with the 
supremacy of Hamilton's principal function. 

In my opinion, there could be no other way to deem a logically consis- 
tent mathematical formalism as inadequate than by demonstrating 
the departure of its consequences from experience or by proving that 
its predictions did not exhaust the possibilities of observation^ pp. 
200-241]. 

We consider this kind of fluxes, because they appear in the classical limit of 
quantum mechanics: 

In the classical approximation, ip describes a fluid of non-interacting 
classical particles of mass m, and subject to the potential V(r): the 
density and current density of this fluid at each point of space are 
respectively equal to the probability density P and the probability 
current j of the quantum particle at that point. 

Indeed, since the continuity equation of this fluid is satisfied [eq. 
(VI. 19)], it suffices to show that the velocity field 

P m 

of this fluid, actually follows the law of motion of the classical fluid 
in question. (SJ p. 223]. 

In this way, this paper demonstrates that there are classical scenarios, as 
the one described by © that cannot be reproduced from quantum mechanics 
in the classical limit. 

In the last section we use our generalization of the Hamilton-Jacobi theory 
to prove — as a mathematical fact, because our proof is logical, not ontological — 
that the wave functions of Klein-Gordon theory can be soundly interpreted as 
describing the motion of ensembles of particles under the action of the electro- 
magnetic field, alone, without any quantum potentials, uninterpreted hidden 
variables, or zero point fields. What Bohm called the quantum potential cor- 
responds in part to kinetic energy and in part to reaction forces associated to 
local creation/annihilation processes. In this form, we address the only objec- 
tion that Einstein raised against the de Broglic-Bohm interpretation of quantum 
mechanics [7] (SJ pp. 33-40]: that, according to Bohm's theory, in the stationary 
states of particles inside a box with rigid walls, the particles had to be at rest. 

The reader should realize that to say that "the wave functions of Klein- 
Gordon theory can be soundly interpreted as describing ensembles of particles 
under the action of the electromagnetic field..." is not the same that to say "the 
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wave functions of Klein- Gordon theory describe ensembles of particles under the 
action of the electromagnetic field..." Therefore, we are not making any onto- 
logical assertions in this paper, which keeps us out of the, otherwise fascinating, 
realm of philosophical speculation. 

We do not say, like Einstein-Podolsky-Rosen, that quantum mechanics is 
not a complete theory of motion|5]. Instead, we say that the set of terms 
that specify the wave function in Madelung's representation, {p, 51201 j can 
be complemented introducing a four vector cj,;, in such manner that the same 
phenomena can be explained using the classical concept of particles that move 
along well defined trajectories. In addition, the logical necessity of introducing 
this four-vector will be mathematically proved beyond any doubt on the basis 
of Hclmholtz theorem. 

This four-vector appears in the equations as the potential of a kind of elec- 
tromagnetic field, that does not produce any Lorentz force on the particles of 
the ensemble under examination. The electric force exactly cancels the mag- 
netic force. However there is not any a priori condition on the divergence of 
the corresponding Faraday tensor. Therefore this field can have an associated 
virtual electric current. 



2 The Helmholtz's Theorem in Three Dimen- 
sions 

To find the best approximation of a vector field / by means of the gradient of 
a potential in a region of space we can follow a variational approach: 
Determine the extreme of the functional: 

F = \j '{V^-ffdV. (6) 

The first variation is: 

SF= f SV4>-(V(t>- f)dV = (7) 
Jn 

I V • [50(V0 - f)¥V - I Sct>V-{Vct>-f)dV. 
Jn Jn 

Making use of Gauss theorem, the first of the last integrals is replaced by 
an integral on Si's boundary(T(S7)). Therefore, 

SF= I 5<j>(y<j) - f)-dS- I 8<PV ■ (V(f> - f)dV (8) 
Jr(n) Jn 

We do not fix the values of <fi on r(S7), so that, from the condition SF = 0, 
we can get a differential equation and a set of boundary conditions: 

V 2 = V ■ /, (9) 
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Vx e r(0) : (V0 - f) ■ h(x) = 0. (10) 

(Where h(x) is the field of unitary normals on r(il).) 

From the theory of harmonic functions we know the last problem has a 
unique solution which we will not discuss further. The field / is thus expressed 
as the sum of two fields: 

f = V<j> + t, 



where t is a solenoidal field, such that 

V x / = V x t, 

and 



(11) 
(12) 



has the minimum value compatible with equation (|12f> . 

Therefore, the field t is a solution for another variational problem: 



5 j n l P+ *-Vx(f-?) 



dV, 



(13) 



where the components of A are Lagrange's multipliers. 

The variation with respect to A leads to the condition V x t = V x /, which 
we already know. 

The variation with respect to t leads to: 



t ■ St - A • V x St 



dV = 0. 



(14) 



Using the identity 

V-(flxl)=l-Vxfl - a-V xb 
to obtain the substitution 

— A • V x St = V • (A x S?) - St- V x A, 
equation (|14f) can be written as: 

(t- V x A) ■ St + V ■ (A x St)] dV = 



(15) 



(t - V X A) • StdV + (n x A) • StdS = 0. 
Jr{fi) 

From this we get the differential equation: 

t = V x A 

From equation (fl2"|) we get: 

V x (V x A) = V x /. 



(16) 
(17) 
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We can choose A so that V • A = 0. For this all we have to do is the 
substitution: 

A -> A + VV>, 

where ip is a particular solution of a Poisson equation 

V 2 ip = -V • A. 

This transformation does not change the fundamental relation l|ltj|) . 
Equation (|17|l is transformed into: 

V 2 A=-Vx/. (18) 

with the extra and boundary conditions: 

V-A = (19) 

Vxe r(f2)(V x X-t) x h = (20) 

Those relations univocally determine the field t. 

We have thus an addition to Helmholtz' theorem in the form: 

Every vector field / with continuous partial derivatives inside a re- 
gion f2 can be written in the form: 

/ = V0 + V x A, 

where A is solenoidal and <j) and A are solutions of the Poisson equa- 
tions: 

V 2 = V • f, 

and 

V 2 A = -V x f. 
If <j> and A are chosen to meet the boundary condition 

Vf G r(fi) : (V0 - /) • h{x) = 0, 

then V</f is the best approximation of the field / as the gradient of 
a scalar function. 

Consider now a vector field ♦••,/»») m a space of n dimensions. As 
before, we find the solutions of a variational problem SF = 0, where 

where 17 is a region of a rt-dimensional space limited by a hyper-surface 
(We are using the summation convention.) 
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The condition SF = leads to a set of corresponding equations, analogous 
to J9j I10|) . As to the vector ti = fi — , we know it is solenoidal, and that it 
is a solution of the variational problem 



S I -UU + K l3 ^—(ti-fi)-— {ti -fi))dV = 0, ( 22 ) 



where is an antisymmetric set of Lagrange's multipliers. 
The extreme condition is easily found to be: 

- Z <»> 

Therefore, U is the divergence of an antisymmetric second order (cartesian) 
tensor. 
Writing 

<S>H = Ay+Jy0, (24) 

we can see that 

meaning that any vector field can be written as the divergence of a second order 
(cartesian) tensor. 

3 Classical Mechanics 

Let Pi (<f 1 , • ■ ■ j 9n, be a vector field, defined in the space-time of a mechanical 
system whose evolution is determined by a Hamilton function 

H(qi,---,q„,pi,---,p„,t), 

so that: 

9H . dH 

% = -5— Pi = —5—- ( 26 ) 
dpi d^i 

(Given that functions defined in Cartan's espace des etats \12\ 174] will be 
projected into functions defined in space-time, to avoid confusions, we will use 
symbols of the kind 

G a p... (27) 
to represent their derivatives. For example: 

dG 



To represent the derivatives with respect to the space-time coordinates we will 
d d_ 

d qi ' dt ■ 



use the usual notation Mr.) 
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We suppose the pi are the momenta of an infinite ensemble of particles, of the 
kind used in Euler's hydro-kinematics. The particle that occupies the position 
qi at time t, gets an increment of action in the interval dt, which is given by: 

dS = ^2 Pi dqi- Hdt. (28) 

The conditions under which this is relation is integrable are 

OQj dqi 

and 



The derivative of pi along the corresponding trajectory in phase space is: 
Pi = -H qi = -^ + ^1i Wj . (31) 



From this we get the equation: 

dpi 



Also, 



Therefore: 



dpi dH_ _ \ - . ,dp L _ £>Pi_. 
dt + dq t ~ ^ qA dq t dqy 



Thus, we get the important conclusion that the vorticity of (pi, ■ ■ ■ ,p n ) is 
zero in the configuration space if and only if the vorticity of (— pi, • • • , — p n , H) 
is zero in space-time. 

According to equation (l.'S2l) 

d 2 pi _ dH qi \ - dH Pk dpi x - <9V; 



s^dHmdPi sr n d2pi (<u\ 



dtdqj dqj ^ dqj dq k ^ k dq k qj 
From this we can see that: 

d f dp, dpA . d ( dpi dp 



dt \dqj dqi J ' ^^ h dq k \dqj dqi) 
dH qj dH qt t y, dH Pk dp j dH Pk dpi 



dq { dqj ^ % dq k dqj dq k 
Considering the relations 

dH <n _ dH i, _ tt dp k jj dp k 

dq t d Ql -^ Hq ^d qi Hq ^dq^ W 
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and 

\p dH Pk dpj dH Pk dpi = 
^ dq, dq k dq 3 dq k 

sr (v 9 Pj _ v ^El\ , rr u ( d P™ dp 3 dp m dp, 

L { M ™ dq k U »« dq k ) + 2^2^ y dq . dqk dq . g Qk 

we can write equation l|35|) in the form 

d f dp, dp A . d ( dpi dp 



dt \dqj dqi J ' ^^ k dq k \dqj %/ ^ 38 ^ 

V tt ( dpj dpk \ + h (f*Eh._ ?Ei\ + 

2^ H ^\n„ k d qJ +Hp ^{dq t dq k l + 



/ dp m dpj dp m dpj 

2^ 2^ p«P™ [ dqi g qk g g qi 



'Jj <Jq k 

Now we assume that the mass matrix is diagonal, so that 



v^v^„ / dp m dpj dp m dpi \ \ - / dp k dpj dp k dpi 

2^2^ Vu Vm \ g q . g qk g q . g q J - 2^ mk [ dq . dqk 9q . dqk 



(39) 

\ - / dp k dpj dpk dpj 
^ k \ dqi dq k dqj dq k 

X^m ( dpk dpi dpi dpj | dpi ®2l dpk dpi 
^ k \ dqi dq k dq k dq k dq k dq k dq 3 dq k 

^ \ dqi dq k dq k dq k dqj dq k/ 

Equation 1(38(1 gives us the material derivative of the vorticity of the field pi. 
From equation 1)89(1 we can see that, under this assumption — that the matrix of 
masses is diagonal — , if the vorticity of (p%, ■ ■ ■ ,p n ) is zero at t = to, it will be 
zero at any other time. Therefore, in those cases, there is a function defined 
in space-time such that: 

d$ TT d<l> 

Obviously, $ is a solution of the Hamilton Jacobi equation: 

9$ <93> d$ 
H( qil ---,q n ,— ,t) + — =Q. (41) 
dqi dq n dt 

Here we are assuming that the trajectories of the particles in the ensem- 
ble do not cross each other, which cannot be granted from the sole condition 
on the vorticity. Therefore we can expect that the solution of (f 4 If) will have 
singularities, depending on the initial conditions. 
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But suppose that {p\, - ■ ■ ,p n ) is not a potential field. The vorticity of 
(—pi, ■ • ■ , — pm H) will not be zero and, according to our previous conclusions, 
there is at least one field (— A±, ■ ■ ■ , —A n , 6) such that: 

and (—px—Ax, • • • , —p n —A n , H+O) is the best approximation of (— J>i, ■ • • , — p n , H) 
as a potential field. Consequently, there is at least one function <!> such that 

pi = ¥- -Ai>B = -5r-e- (43) 

dq. t dt 
The function $ is a solution of the differential equation 

H(qx, ■ ■ ■ , q n , ^ Ai, • • • , A„, t) + 6 + — = 0, (44) 

aqi dq n at 

which is the equation we had obtained, had we started with a potential field 
(j>x, ■ ■ ■ ,p n ) and another Hamilton's function: 

H'(qx, ■ ■ ■ ,q n ,px, ■ • ■ ,Pn,t) = H(qx, ■ ■ ■ ,q n ,Pi - M, ■ • • ,Pn - A n ,t) + 9 (45) 

It is clear that equation l|44l) is not enough to determine the fields $ and A, 
but this difficulty is easily surmounted. To show this in a way that is free of 
mathematical complexities, we will consider the case of a single particle, where 
the Hamilton's function of the original system is: 

-Q 

H(p,q) = ^-+V. (46) 
The new Hamilton's function is 

H'(q,p) = ^^ + V + Q, (47) 
2m 

which is formally analogous to the Hamilton's function of a particle in an elec- 
tromagnetic field. The corresponding Hamilton-Jacobi equation is: 

9$ 

H'(q,V$) + — =0. (48) 

From l|47ll we see that the function and the vector A have to be chosen so 
that the corresponding Lorentz force is equal to zero for the given field: 

-V6-^ + ^(£W)x(Vxl)=0. (49) 

Equations l|48() and (|49|1 constitute the generalization of Hamilton-Jacobi 
theory to the consideration of fields with vorticity. 
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To make an example, we have considered an ensemble of free particles that 
satisfy the initial conditions: 

v(fo,t ) = u> x ro, (50) 

where to = Lok is a constant vector in the direction of the z axis. The position 
at time t of the particle that occupies position ro at time to is 

f(t) =r + {ujxr )(t-t Q ). (51) 

Or, in open form: 

x(t) \ I 1 -w(i-t„) \ / x 
y(t) = [ Lj(t-t ) 1 T/o | . (52) 

z(t) J V 1 / V z o 

The determinant of this system is greater than zero. Therefore, the transfor- 
mation is invertible. Actually, we have this: 

x + uj(t-t Q )y 

, .s y-us(t-t )x 

Vo{*,V,*,t)= 1 + b p {t _ tor 

z (x,y,z,t) = z. 
From this and l|51|) we can get the field of momenta at time t: 

-mwy + muj 2 (t-t )x 
Px = -mojyo = 1 + Lu2(t _ to)2 (54) 

mux + muj 2 (t — to)y 

Py = mujXo = i+^(t- to r ' 

p* = o, 

and the field of kinetic energy 

1 mw 2 {x 2 + y 2 ) 

21 +uj2(t- to) 2 ' ( j 

The potential and solenoidal parts of the field of momenta are clearly sepa- 
rated. 

p = V<P-A, (56) 

where 

lmuj 2 (t-t )(x 2 +y 2 ) 

2 i + Lu 2 (t-t y 



± hum yu - uojyu, T y 



Ax ~ i+u*(t-t r (58) 
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mux 

~ y = l+uj^t-to) 2 ' 

and 

-A z = 0. 

Also: 

dt l + u 2 {t~t ) 2 y ' 

At this moment we can see that the alternative Hamilton's function for this 
problem is: 



2m 



mojy \ I mux \ 2 2 

Pv+ 1 , , ar* *.\» + P» 



■ l + w 2 (i-io)V V V l + ^(t-t Q )\ 

(60) 

muj 2 (x 2 + y 2 ) 
~l+u 2 {t-t ) 2 ' 

Because of the way we have constructed this function, it is not difficult to 
see that <|57[) is a particular solution of the Hamilton- Jacobi Equation 



2/ \2/\2 

d$ mujy \ / <9<& muix \ 1 <9<I> " 



(61) 



2m l. U.r I f ^-(/ - /,,)-' y ' V /V// I + / - /„)- 

mw 2 (a; 2 + w 2 ) 3$ 
H = 

i + co 2 (t-f ) 2 at 

That H49|) is also satisfied can be easily proved from the equalities: 

^ n 9A _ 2muj 2 (xi + yj) 2mu 3 (t - t )(yi - xj) 

dt ~ l+u 2 (t-t a ) 2 (i +w 2^_ to)2 )2 > I > 



mu) 2 (t — to)x — toy * | lux + muj 2 {t — t[)) 2 y 



l + Lj 2 (t-t ) 2 ) " ' l + ^ 2 (t-t ) 2 ) j (63) 
and 

^ ? —2mujk , 

VxA= — . (64) 

l + cj 2 (<-t ) 2 v y 

Passing to special relativity, in previous papers and 0, we have observed 
that in the case of a field of four-velocities that represents the motion of an 
infinite ensemble of particles, the derivative of the four-velocity along the cor- 
responding world-lines 

duj „ du; 



ds dxi 



(65) 



is determined by its vorticity. This is so because, from the condition vPiij = 1, 
it follows that 

i du 3 A 
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Combining this and i|t)5|) we get 

d7 = u W-W' (66) 

Then — we see now — the vorticity of a field of four-velocities cannot be zero 
at least it describes an ensemble of free particles. If the particles move under 
the action of an electromagnetic field we have: 

where 

P L = mcu t Ai, (68) 

c 

m is the rest mass of the particles, and q is the corresponding electric charge. 
We can write the vector Pi in the form 

9S 

Pl = ~d^ + ^ (69) 

From JflSJ), lESt ) EOt : CP ) we get a non-linear system of differential equa- 
tions: 

9S . 1 4k , ,.k\ f 9S , 1 A , , , \ _ ™2 2 



where, because of (IHTIl : 



and 



' c 1A * + ^l-^ + ^ + <*)="> v ' < 72 > 

These equations encompass the most reasonable generalization of Hamilton- 
Jacobi theory to special relativity since, as we said before, a potential field of 
four-velocities represents necessarily a field of free particles. Vorticity plays a 
special role in relativistic Hamilton-Jacobi theory, and there is not a physical 
reason to believe that the only real fields of four- velocities are those for which 
there is a scalar function (j) such that ^2 p. 488-509] 

mcv u_g A u = d^ (?5) 
C OX u 

Quite the contrary, this reformulation of Hamilton-Jacobi theory allows us to 
prove that it is possible to interpret the wave functions of Klein-Gordon theory 
as describing the motion of an ensemble of particles under the action of the elec- 
tromagnetic field, alone, without quantum potentials or uninterpreted hidden 
variables, where the number of particles is not locally conserved. 
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4 The Klein-Gordon Field 



From the Klein-Gordon equation 

d 2 , 2ih ... d^> 



-h' - — * A"— + 9 i #4 /1 * = mV$, (76) 

ax^ox^ c oxf* 

we can easily show that 

d (%h, „ <9*\ g , . \ 

w— l-s-C* ) - = 0. 77 

dx M V 2 ^ ^ c / 

From this, using the Madelung substitution 

* = ^pe {t/h)s , (78) 

we prove that 

This suggest that p could be interpreted as a density of particles, in the 
system of reference where they are at rest and that 



v» = - 



1 dS q 



mc dx^ mc 



1^, (80) 



which is usually rejected on the grounds that is not unitary time-like by 
definition. From the previous section of this paper we see now that we can 
complete the representation if we suppose that 

dS q . . 



where 



and 



The continuity equation l|79|) is transformed into: 

dx^ dx^ ' 



(84) 



that shows that the vector w M as linked to local production/anihilation processes. 
This is not too strange if we consider that, after all, this four-vector represents 
a kind of electromagnetic field. Equation Ij82(l is the condition it does not 
produce a Lorentz' force. However, its very existence, implies the appearance 
of the corresponding conserved current. 
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From equations 17(i|) and (|81|l we can show that 



-2mc,<H + u^ = -^VP (85) 

Equations Ij81|l to l|85|l make up a complete description of an ensemble of 
particles that move under the action of the electromagnetic field, in such way 
that the number of particles is not locally conserved. Equations Ij84(l and H85[l 
are the only ones that include the density of particles and Planck's constant Ti, 
which is then interpreted as an empirical parameter that determines the local 
rate of creation/anihilation of matter. The reinterpretation of Dirac's field is 
similar to the one we have exposed in 0. 

In the low speed limit, we replace the four- vector uji by its best approxima- 
tion in the sense we studied in the second section. This approximation does 
not meet the requirements of Lorentz invariance, but Lorentz invariance is not 
required in the low speed limit. Then we can follow the same line of though we 
followed in pQ and to recover Schrodinger's and Pauli's equations. 
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